The potential energy of a 40 K Fermi gas in the BCS-BEC crossover 
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We present a measurement of the potential energy of an ultracold trapped gas of 40 K atoms in the 
BCS-BEC crossover and investigate the temperature dependence of this energy at a wide Feshbach 
resonance, where the gas is in the unitarity limit. In particular, we study the ratio of the potential 
energy in the region of the unitarity limit to that of a non-interacting gas, and in the T — limit we 
extract the universal many-body parameter /3. We find /3 = — 0.54^q'^2; this value is consistent with 
previous measurements using 6 Li atoms and also with recent theory and Monte Carlo calculations. 
This result demonstrates the universality of ultracold Fermi gases in the strongly interacting regime. 
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With the emergence of novel Fermi gas systems, ex- 
perimentalists can now access the Bardeen Cooper Schri- 
effer (BCS)- Bose-Einstein Condensate (BEC) crossover 
in ultracold gases of atoms. Using atomic scattering 
resonances in gases of 40 K and 6 Li, it is possible to 
widely tune the interatomic scattering length, a, and 
move continuously between a gas of weakly interact- 
ing fermions and a gas of condensed molecules. The 
BCS-BEC crossover occurs in the strongly interacting 
regime where the scattering length is large enough that 
— 1 < 1/kpa < 1, where kp is the Fermi wave vector. 
Experimental studies of these strongly interacting Fermi 
systems have revealed many interesting properties includ- 
inga phase transition involving condensates of atom pairs 
[3j3] , a pairing gap [3] , and vortices in a rotating gas [1] . 

As a gas of fermions is cooled from the classical regime 
to quantum degeneracy, the Pauli exclusion principle be- 
comes manifest in the properties of the ultracold gas 
For example, a zero-temperature Fermi gas in 
a confining potential has a finite energy and a finite size 
due to Fermi pressure, which is responsible for the sta- 
bility of white dwarf and neutron stars. As the two-body 
scattering length is tuned to be arbitrarily attractive, one 
would expect that the gas should be compressed due to 
attractive interactions and pairing effects 0, H| • 

This behavior has received theoretical consideration 
and should not depend on the details of the inter- 
atomic potential for a wide Feshbach resonance, where 
the Fermi energy is much smaller than the energy equiv- 
alent width of the resonance [3, [13, [HI, [e| • Furthermore, 
at resonance, where the two-body scattering length a di- 
verges, the energy of the gas is expected to be universal 
[l3L [l4 [1^ in that it depends only on the Fermi en- 
ergy Ep and the relative temperature T/Tp. The den- 
sity profile of a trapped T — unitarity limited gas is 
then expected to be simply a rescaled version of the non- 
interacting density profile. This results in a simple rescal- 
ing of the size and energy, which can be parameterized 
by a universal many-body parameter j3 3, 13] • 

Experimentally, (5 has only been reported for 6 Li where 
the size and e nergy of a trapped gas has been examined 
@, US El E3, Hfl HE S3- The most precise determina- 



tion, fi = — 0.54±0.05, was reported recently in Ref. [21 1. 
While this value is in good agreement with the predicted 
value, a measurement using a different atomic species 
is essential to demonstrate the universality of strongly 
interacting Fermi gases. Here we report on a measure- 
ment of the potential energy at resonance for an ultracold 
gas of 40 K. We have also measured the potential energy 
throughout the strongly interacting regime and investi- 
gated the temperature dependence of the potential en- 
ergy at resonance. 

For these experiments we cool a gas of fermionic K 
atoms to ultracold temperatures using previously de- 
scribed methods (23|. A nearly equal mixture of the two 
lowest energy hyperfine spin states, ]/, m/) = 1 9/2 , — 9/2} 
and |9/2,— 7/2), is confined in a crossed beam optical 
dipole trap. The trap consists of a horizontal laser beam 
parallel to z with a \ radius of 32 /j,m and a vertical 
beam parallel to y with a j% radius of 200 /im. For the 
experiments reported here the harmonic trap frequencies 
were typically uv/27r = 184 Hz and uj z /2it = 18 Hz. Ap- 
proximately 10 5 atoms per spin state are cooled to a final 
temperature of T ps 0.08 Tp, where Tp = Ep/kf, is the 
Fermi temperature and kb is Boltzmann's constant. 

The final evaporative cooling of the gas occurs on the 
BCS side of the magnetic-field Feshbach resonance where 
a « —1000 clq. The optical trap is then ramped up to ap- 
proximately 1.5 times of the shallowest trap depth used 
for evaporation. To vary a, the magnetic field is adiabat- 
ically ramped to various final values. The optical trap 
is then suddenly switched off and the gas is allowed to 
expand for 1.867 ms. During this short expansion time 
there is significant expansion in the radial direction but 
negligible expansion of the cloud in the axial direction. 
We then use absorption imaging to probe the density dis- 
tribution of the atom cloud. The probe beam propagates 
along one of the radial directions, x, and is pulsed on 
for 40 jis. For each absorption image we perform a 2D 
surface fit to a finite temperature Fermi-Dirac function 
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where £, a y , a z , and pk are independent fitting parame- 
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FIG. 1: Measured potential energy E pot normalized to the 
value measured in the non-interacting regime Ep Bt versus 
1/kpa. Here (°) denotes a quantity measured in the non- 
interacting regime, i.e., at the zero crossing of the s-wave 
scattering length. The resonance is located at 202.10 ± 0.07 
G the dashed lines show the uncertainty in the resonance 
location. Data points toward the BCS limit show good agree- 
ment with a zero temperature mean-field calculation (solid 
line) . The larger error bars on the BEC side of the resonance 
reflect uncertainties due to heating of the gas due to inelas- 
tic loss. In the strongly interacting region their exists ±0.1 
uncertainty in l/k° F a due to uncertainty in the resonance po- 
sition, (inset) Subset of the data focusing on the strongly 
interacting region near resonance. 



ters and <7 n (x) = pr UM- This is the expected optical 
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depth (OD) distribution for a non-interacting cloud both 
in trap and after expansion. Empirically we find that this 
function also fits well in the strongly interacting regime. 
The potential energy of the trapped gas is obtained from 
the cloud profile in the axial direction. The potential 
energy per particle in the axial direction (z) is given by 
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where m is the mass of 40 K. 

It is useful to normalize the measured potential energy 
of the strongly interacting gas to that of an ideal (non- 
interacting) Fermi gas. In previous experiments using 6 Li 
atoms, the measured cloud sizes and energies were nor- 
malized to a calculated value for the non-interacting gas. 
This can introduce systematic errors because the calcu- 
lation relies on the atom number and trap frequencies, 
which can have systematic errors. However, for the 40 K 
Feshbach resonance, we are able to reduce the system- 
atic uncertainty by measuring the potential energy of the 
non-interacting gas. This can be accomplished by going 
to the zero crossing of the resonance where a — 0, which 



FIG. 2: Measured potential energy at the Feshbach resonance 
versus magnetic-field ramp duration. For very fast ramps, we 
measure a higher energy because of nonadiabaticity. For very 
slow ramps, heating due to inelastic collisions increases the 
measured energy. 



is only 10 G away. The measured potential energy ratio 
for an ultracold 40 K gas is displayed in Fig. Q]as a func- 
tion of 1/kpa. In this paper we use a superscript naught 
( ) to denote measurements made in the non-interacting 
regime. The inset to Fig. [T] focuses on the strongly inter- 
acting region near the resonance. As expected, the data 
show that the interactions cause a strong reduction in the 
potential energy due to a compression of the trapped gas. 
One would expect that on the BEC side of the resonance 
Epot will depend on condensate fraction. For tempera- 
tures similar to these experiments we find a maximum 
condensate fraction of approximately 15% on resonance; 
this fraction decreases as detuning from the resonance is 
increased [l[. 

The error bars in Fig. [1] include statistical uncer- 
tainty in repeated measurements as well as an uncer- 
tainty due to heating during the magnetic-field ramps. 
The magnetic-field ramps must be sufficiently slow to be 
adiabatic; however, heating during the ramp can be a 
problem for slower ramps. For the different final mag- 
netic fields we investigated the dependence of the mea- 
sured potential energy on the duration of the linear ramp. 
An example of this is shown in Fig. [21 where the gas 
was ramped from the magnetic field used for evaporation 
(203.39 ± 0.01 G) to the resonance position. 

To determine the optimum ramp rate, as well as the 
effect of heating on the potential energy measurement, we 
fit data such as that shown in Fig. [5] to an exponential 
decay plus linear heating. From the fit we determine 
the final potential energy of the cloud if heating were 
not present. This introduces a correction that is applied 
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to the data shown in Fig. [TJ Note that on the BCS 
side of the resonance we see little or no heating due to 
magnetic-field ramps, and the error bars are dominated 
by shot-to-shot statistical uncertainty. 

We can gain some theoretical insight into the effect of 
interactions on the energy of our trapped gas by consid- 
ering a simple mean-field approach. While this approach 
neglects pairing and therefore is not sufficient to fully un- 
derstand the behavior of our gas, it provides a flavor of 
how the potential energy is affected by interactions. Fol- 
lowing the argument outlined in Ref. [rjj, the equation 
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of state for an ideal zero-temperature Fermi gas is 

fi = e F (x) + U MF {x) + U trap (x), 

where [i is the chemical potential, £f(x) is the local 
Fermi energy, £/mf(x) is the mean- field contribution, 
and C/( rap (x) is the trapping potential. We can relate 
e F (x) = 7^k F (x) to the density n(x) = g^fc|(x) via 

£f(x) = ^■[67r 2 n(x)] 2 / 3 . The interactions appear in the 
density dependent mean-field contribution, [/mf(x) = 
^ 2 ^- 2 -n(x). This equation can be solved self-consistently 
to determine the in-trap density profile of the cloud, and 
thus the potential energy per atom. In Fig. [TJ we com- 
pare the data to the mean-field calculation (solid line ) fo r 
the normal state on the BCS-side of the resonance [231 ] . 
Near the resonance, in the strongly interacting regime, it 
is clear that this approximation breaks down and a more 
sophisticated theory is required. 

Very near the resonance the scattering length a di- 
verges and the equation given above for [/mf(x) becomes 
unphysical. At resonance (l/k F a = 0), the only energy 
scale is the Fermi energy; this give us an approximate 
effective scattering length a e ff = — \/k F . This substi- 
tution shows that the local mean-field energy is propor- 
tional to the Fermi energy, and one can define a constant 
of proportionality [3 given by C/mf(x) = f3e F (-x) [si Il6j|. 
In this simple mean-field estimate, Umf{*) = — ^cf(x), 
or /3mf = —0.41. The negative sign for the scattering 
length is not obvious from this approach, but a more 
sophisticated many-body approach shows the mean-field 
interaction should be attractive 
Eq. © as 



131 ]. Now we can write 



H= (l + f3)e F (x) + U trap {x). 



(4) 



Solving for the density profile for a harmonic trap and 
then integrating to find the energy per particle, one finds 
that the potential energy of a T = gas in the unitarity 
limit is simply E pot = just as in the case of a non- 
interacting Fermi gas. To find the ratio of the chemical 
potential at resonance fi to that for a non-interacting gas 
fjP, we hold the number constant for each case, N 



• a (x)d 3 x = J n(x)d 3 x, to find 



UJ 



1.1 



1.0 



* 0.9- 



0.8 
0.7 
0.6 




/<>' 



(5) 



FIG. 3: Potential energy E pot normalized to the measured en- 
ergy in the non-interacting regime Ep 0t vs the non-interacting 
gas temperature (T/T F )°- The cloud is heated by parametri- 
cally modulating the trapping potential. For these data the 
trapping frequencies in the radial direction vary from ~180 
Hz to 450 Hz and in the axial direction from ~18 Hz to 21 
Hz. 



Thus, the universal parameter (3 can be extracted by 
measuring the ratio of the potential energy at resonance 
to the potential energy of a non-interacting, trapped 
Fermi gas, in the T — limit. 

From the data in Fig. [TJ we find E/E = 0.70 ± 0.02 
on resonance giving (3* = —0.51 ± 0.03. Whereas (3 
is normally only defined for T = 0, we introduce (*) 
to denote that the system is at a finite temperature 
{T/T F f = 0.08 ± 0.01. This experiment was conducted 
at the Feshbach resonance and the error bars include 
statistical error as well as the heating effects mentioned 
above. Including uncertainty in the resonance position 
we find (3* = -0.51±°;5*. 

We now consider experimentally the question of 
whether the gas is sufficiently cold to be in the T = 
limit as required for an accurate determination of (3. In 
Fig. [3] we show the measured potential energy ratio as a 
function of (T/T F )°. In this experiment, we heat the gas 
by parametrically modulating the optical trap strength. 
However, before heating the gas, we first increase the 
optical trap depth to prevent number loss and ensure 
harmonic confinement. Both the resulting temperature 
(T/T F )° and k F are determined using an ultracold cloud 
prepared as described above and slowly ramping the mag- 
netic field to the point where a — just before the trap 
is turned off. The gas is allowed to expand freely for 14 
ms and fit according to Eq. |TJ. We extract {T/T F )° 
from the fugacity £ of the non-interacting cloud using 
92(— C) = — (T/T F )~ 3 /6. It is important to note that in 
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this experiment we expect the magnetic-field ramps to 
keep the entropy constant but not the temperature. 

The data in Fig. [3] clearly show that the universal 
many-body parameter (3 depends strongly on tempera- 
ture. Furthermore, it is not clear that a gas at (T/Tp)° — 
0.08 is sufficiently cold to determine the T = limit of 
(3. For the purpose of extrapolating to zero tempera- 
ture, we fit a quadratic function to the data points below 
(T/T»° = 0.25 for which we find [3 = — 0.54±g;?|. The 
error bars reflect the uncertainty in the extrapolation to 
T = and the uncertainty in the resonance position. 
This value of the universal many-body parameter (3, as 
well as the value at (T/Tp)° = 0.08 .is in good agree- 
ment with Monte Carlo calculations 25l 26, 

" 31 
30]. 



m par- 
These values 



and recent theoretical calculations 30, 
ticular, we note (3 — —0.545 in Ref. 
are also in good agreement with multiple experimental 
reports in 6 Li: (3 — 




0.61 ±0.15, -0.49, and 



0.54 ± 0.05 in Refs. [3|, [20j, |22j, and |2l| respec- 
tively. We also note that from the kinetic energy mea- 
surement in 40 K reported in Refs. 23, 35|, we can extract 
f3 = —0.62 ± 0.07 j36|. This is in good agreement with 
the value of (3 found using the potential energy presented 
in this paper. 

In summary we have studied the potential energy of 
a strongly interacting quantum degenerate gas of K 
Fermi atoms. At resonance limit our results are con- 
sistent with current theory as well as previous experi- 
ments in 6 Li, thereby strengthening the theory of univer- 
sality of these Fermi gas systems. Universality necessar- 
ily assumes that the s-wave Feshbach resonances in 40 K 
and 6 Li are wide. The question of whether or not the 
40 K Feshbach resonance is wide has been under debate 
lol . 12l 37 . 38j | : however, the good agreement of our re- 
sults with 6 Li provides compelling evidence that 40 K is 
also a wide resonance. We have also measured the tem- 
perature dependence of a universal many-body parame- 
ter that could be compared to recent Monte Carlo results 
for the temperature dependent energy of a homogeneous 
Fermi gas j29ll39jj. 
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